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Abstract 

Improving our understanding of the initial conditions 
and earliest stages of protostellar collapse is crucial to 
gain insight into the origin of stellar masses, multiple 
systems, and protoplanetary disks. Observationally, there 
are two complementary approaches to this problem: (1) 
studying the structure and kinematics of prestellar cores 
observed prior to protostar formation, and (2) studying 
the structure of young (e.g. Class 0) accreting protostars 
observed soon after point mass formation. We discuss re- 
cent advances made in this area thanks to (sub)millimeter 
mapping observations with large single-dish telescopes 
and interferometers. In particular, we argue that the be- 
ginning of protostellar collapse is much more violent in 
cluster-forming clouds than in regions of distributed star 
formation. Major breakthroughs are expected in this field 
from future large submillimeter instruments such as Her- 
schel and ALMA. 
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1. Introduction 

Although the formation of low-mass stars is now rea- 
sonably well understood in outline (see, e.g.. Mannings, 
Boss, & Russell 2000 for recent reviews), several impor- 
tant aspects remain poorly known, such as the initial 
stages of the collapse process, the mechanism(s) select- 
ing stellar masses, or the formation of multiple systems. 
Some progress on the earliest stages of star formation was 
achieved over the last decade thanks to the use of sensitive 
receivers on large (sub)millimeter radiotelescopes such as 
JCMT, CSO, and the IRAM 30m. Young protostars were 
identified at the beginning of the main accretion phase 
(Class objects - Andre, Ward-Thompson, & Barsony 
1993), and the starless dense cores extensively studied 
in NH3 by Myers and collaborators (e.g. Benson & Myers 



1989) were found to be characterized by flat inner density 
gradients (Ward-Thompson et al. 1994). Direct evidence 
for infall motions was observed toward a large number of 
Class protostars and prestellar cores (e.g. Gregersen 
et al. 1997 - see § 3 below and Evans, this volume). 
Class objects were also found to drive more powerful 
outflows than more evolved (Class I) protostars, suggest- 
ing a marked decrease of the mass accretion/ejection rates 
in the course of protostellar evolution (Bontemps et al. 
1996). As advocated by Henriksen, Andre, & Bontemps 
(1997), there may be a causal relationship between these 
results, in the sense that the accretion/ejection decline 
during the protostellar phase may be a direct consequence 
of the form of the density profile at the prestellar stage. 
Further studying the detailed properties of prestellar cores 
and young protostars is of prime importance to distin- 
guish between collapse models and shed light on the ori- 
gin of stellar masses. Indeed, the effective reservoirs of 
mass available for the formation of individual stars may 
be largely determined at the prestellar stage, and it is 
during the main protostellar accretion phase that stars 
accrete some fraction of these reservoirs and build up the 
masses they will have on the zero-age main sequence. 

After a brief introduction of the theoretical background 
(§ 1.1), we review recent observational advances concern- 
ing the density and velocity structure of cloud cores in 
§ 2 and § 3, respectively. We conclude in § 4 with a com- 
parison between observations and theoretical models. 

1.1. Collapse Initial Conditions: Theory 

The inside-out collapse model of Shu (1977), starting 
from a singular isothermal sphere (SIS) or toroid (cf. Li 
& Shu 1996, 1997), is well known and underlies the 'stan- 
dard' picture of isolated, low-mass star formation (e.g. 
Shu, Adams, & Lizano 1987). 

Other collapse models exist, however, which adopt differ- 
ent initial conditions. In particular, Whitworth & Sum- 
mers (1985) have shown that there is a two-parameter 
continuum of similarity solutions to the problem of isother- 
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Figure 1. Dust continuum maps of L1544 (o-) o,'>^d 
IRAM 04191 (h) at 1.3 mm taken with the IRAM 30 m 
telescope and the MPIfR bolometer array (from Ward- 
Thompson et al. 1999 and Andre et al. 1999, respec- 
tively). Effective resolution: 13" ; base contour and con- 
tour step: 20 mjy/beam. The direction of the magnetic 
field measured in L1544 with the SCUBA polarimeter on 
JCMT (Ward-Thompson et al. 2000) and the collimated 
C0(2-l) bipolar flow emanating from IRAM 04191 are 
shown. 



mal spherical collapse. One of the parameters measures 
how close to hydrostatic equilibrium the system is ini- 
tially, while the other parameter reflects how important 
external compression is in initiating the collapse. In this 
continuum, the solutions proposed by Shu (1977) and 
Larson {1969)-Penston (1969) represent two extreme lim- 
its. All of the similarity solutions share a universal evo- 
lutionary pattern. At early times (i < 0), a compression 
wave (initiated by, e.g., an external disturbance) prop- 
agates inward leaving behind it a p(r) oc density 
profile. At t = 0, the compression wave reaches the cen- 



ter and a point mass forms which subsequently grows by 
accretion. At later times {t > 0), this wave is reflected 
into a rarefaction or expansion wave, propagating out- 
ward through the infalling gas (at the isothermal sound 
speed as), and leaving behind it a free-fall p{r) oc r~^-^ 
density distribution. Several well-known features of the 
Shu model (such as the expansion wave) are thus in fact 
common to all solutions. The various solutions can be 
distinguished by the absolute values of the density and 
velocity at i ~ 0. In particular, the Shu (1977) solution 
has p{r) = {a,^/2T:G) and is static (v — 0) a,t t — 0, 
while the Larson-Penston (1969) solution is ^ 4.4 times 
denser and far from equilibrium {v ~ —3.3 aj. During the 
accretion phase (t > 0), the infall envelope is a factor ^ 7 
denser in the Larson-Penston solution. Accordingly, the 
mass infall rate is also much larger in the Larson-Penston 



case (~ 47 a,^/G) than in the Shu case 



VG). 



In practice, however, protostellar collapse is unlikely 
to be strictly self-similar, and the above similarity solu- 
tions can only be taken as plausible asymptotes. More 
realistic initial conditions than the SIS are provided by 
the so-called 'Bonnor-Ebcrt' spheres (e.g. Bonnor 1956), 
which represent the equilibrium states for self-gravitating 
isothermal spheres and have a flat density profile in their 
central ^ Jeans length region. Such spheres are stable for 
a center-to-edge density contrast < 14.3 and unstable for 
a density contrast > 14.3 (e.g. Bonnor 1956). Numeri- 
cal hydrodynamic simulations of cloud collapse starting 
from such initial conditions (e.g. Foster & Chevalier 1993, 
Hennebelle et al. 2002) find that the Larson-Penston sim- 
ilarity solution is generally a good approximation near 
point-mass formation {t = 0) at small radii, but that 
the Shu solution is more adequate at intermediate t > 
times, before the expansion wave reaches the edge of the 
initial, pre-coUapse dense core. In general, the mass ac- 
cretion rate is thus expected to be time-dependent. 

Observationally, it is by comparing the (density and 
velocity) structure of prestellar cores such as L1544 (see 
Fig. ||a) with the structure of the envelopes surround- 
ing Class protostars such as IRAM 04191 (cf. Fig. Qd) 
that one may hope to constrain the initial conditions for 
collapse and to discriminate between the various existing 
models. 

2. Density Structure 

2.1. Prestellar cores 

Two main approaches have been used to trace the density 
structure of cloud cores: (1) mapping the optically thin 
(sub)millimeter continuum emission from the cold dust 
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contained in the cores, and (2) mapping the same cold 
core dust in absorption against the background infrared 
emission (originating from warm cloud dust or remote 
stars) . Mapping the molecular gas component is generally 
less effective as most molecules tend to freeze out onto 
dust grains in the dense, cold inner parts of cloud cores 
(e.g. Kramer et al. 1999, Walmsley et al. 2001, Bacmann 
et al. 2002, and Bergin 2002). 

Ward-Thompson et al. (1994, 1999) and Andre et al. 
(1996) employed the first approach to probe the struc- 
ture of prestellar cores (see also Shirley et al. 2000). Un- 
der the simplifying assumption of spatially uniform dust 
(temperature and emissivity) properties, they concluded 
that the radial density profiles of isolated prestellar cores 
were not consistent with the single p{r) oc power law 
of the SIS but were flatter than p(r) oc r^^ in their inner 
regions (for r < Rfiat), and approached p{r) cx r^^ only 

2500-5000 AU. The de- 



beyond a typical radius Rjiat 
tails of this conclusion have been challenged by ttSvans et 



al. ( POOI)I and ^ucconi et al. (2001)| who made the impor- 
tant point that starless cores are not strictly isothermal. 
Accounting for a realistic temperature distribution inside 
the cores (see § 2.3 below), these authors found smaller 
values for Rfiat and even claimed that in some cases the 
observed profiles could not be distinguished from a SIS 
profile. We will return to this point below (e.g. § 2.3). 

More recently, it has been possible to use the absorp- 
tion approach, both in the mid-IR from space (e.g. Bac- 
mann et al. 2000, Siebenmorgen & Kriigel 2000) and in 
the near-IR from the ground (e.g. Alves et al. 2001). In 
particular, using the ISOCAM infrared camera aboard 
the ISO satellite, Bacmann et al. (2000) carried out a 
7 /zm survey of 24 low-mass starless cores (all undetected 
by IRAS) and observed absorption features in 23 of them. 
As an example. Fig. |^ shows the ISOCAM image ob- 
tained at 6.75 /im for the prestellar core L1689B in the 
Ophiuchus complex (c? ^ 150 pc). The core is seen in 
absorption against the diffuse mid-IR background aris- 
ing from the rear side of the parent molecular cloud (e.g. 
Bernard et al. 1993, see also § 2.3 below). 

The column density structure of the absorbing cores 
can be derived from a simple modeling of the ISOCAM 
observations. The mid-IR intensity measured at projected 
radius r from core center may be expressed as: 

I{r)=hack-e-^^'^ +Ifore, (1) 

where r is the core optical depth (directly related to the 
column density N//^ via the dust opacity) and hack and 
Ifore are the background and foreground intensities, re- 
spectively. Both Iback and Ifore are uncertain but their 
values are well constrained when independent millimeter 
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Figure 2. (a) ISOCAM 6.75 iim absorption image of 
the prestellar core L1689B (see color scale on the right 
in MJy/sr). The 1.3 mm continuum emission map ob- 
tained by Andr et al. (1996) at the IRAM 30 m tele- 
scope is superposed as green contours (levels: 10, 30, 
50 mJy/13" beam), (b) Column density profile of L1689B 
(crosses) derived from the absorption map shown in (a) 
by averaging the intensity over elliptical annuli for a J^CF 
sector in the southern part of the core. The dashed curves 
show the most extreme profiles compatible with the data 
given the uncertainties affecting the absorption analysis. 
The blue, solid curve represents the best fit of a Bonnor- 
Ebert sphere model (embedded in a medium of uniform 
column density), obtained with the following parameters: 
Pc/ Pout — 40 ± 15 (i.e., well into the unstable regime), 
Teff = 50 ± 20 K, Pext/kB = 5 ± 3 X lO^Kcm-^. For 
comparison, the dotted line shows the N oc profile 
of a SIS at T = 10 K. (Adapted from Bacmann et al. 
2000.) 
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measurements of tlie H2 column density in the centers 
and outer parts of the cores are available (see Bacmann 
et al. 2000 for details). 

Using these constraints on Iback and //ore , a range of 
possible column density profiles can be derived for each 
core imaged in absorption with ISOCAM. This is shown 
in Fig. I^b for L1689B, where the crosses display our 'best' 
estimate of the column density profile (derived using a re- 
vised estimate of 8 x 10^^ cm^^ for the central column 
density from millimeter continuum data, in agreement 
with a central temperature of ^ 8 K - see § 2.3 below - 
instead of 12.5 K as assumed by Bacmann et al. 2000). 
The dashed curves lying above and below the 'best' pro- 
file show the most extreme column density profiles that 
are consistent with the mid-IR absorption data, given the 
range of permitted values for Iback and I fore- In general, 
four different regimes can be distinguished on the prestel- 
lar column density profiles (see, e.g.. Fig. ||b): 

a) a flat inner region (of radius Rfiat = 5000±1000 AU 
for L1689B according to the fltting analysis shown in 
Fig. ^), b) a region roughly consistent with A^h2 oc f^^ 
(corresponding to p oc for a spheroidal core), c) an 
edge where the column density falls off typically more 
rapidly than A^h2 oc with projected radius (suggest- 
ing a density gradient steeper than p oc r~'^), until d) the 
end of the core is reached (at Rout = 28000 ± 1000 AU 
for L1689B) and A^h2 fluctuates about the mean value 
Nout ~ 1 — 2 X 10^^ cm~^ characterizing the ambient 
molecular cloud. 

A qualitatively similar column density profile was de- 
rived by Alves et al. (2001) for the Bok globule B68 (for 
which Rfiat « 4000 AU and Rout ~ 12500 AU), based 
on VLT and NTT measurements of the near-IR colors of 
background stars seen through, and suffering extinction 
from, this globule. 

Thus, the inner flattening of the density profiles of 
prestellar cores first observed in submillimeter emission 
(see above) is confirmed by infrared absorption/extinction 
measurements, which are essentially independent of the 
temperature distribution. Moreover, a new result emerges: 
at least in some cases, such as L1689B and B68, isolated 
prestellar cores appear to be characterized by sharp edges 
defining typical outer radii Rout ^ 20000 AU ^ 0.1 pc. 

2.2. Comparison with theoretical models 

These observed features, i.e., fiat inner region and sharp 
outer edge, set constraints on possible models of core 
structure. First of all, it is clear that self-similar, singu- 
lar models such as the SIS or magnetized, non-spherical 



generalizations of it (e.g. Shu 1977, Li & Shu 1996), can- 
not account for the detailed density structure of prestellar 
cores. As pointed out by Bacmann et al. (2000), the lo- 
gotropic models of McLaughlin & Pudritz (1996), which 
behave as p oc r^^ over a wide range of radii, are also 
inconsistent with the observations of (low-mass) starless 
cores. On the other hand, the circularly-averaged col- 
umn density profiles can often be fitted remarkably well 
with models of pressure-bounded Bonnor-Ebert spheres, 
as first demonstrated by Alves et al. (2001) for B68. This 
is also the case of L1689B as illustrated in Fig. ||b. The 
quality of such fits shows that Bonnor-Ebert spheres pro- 
vide a good, first order model for the structure of isolated 
prestellar cores. In detail, however, there are several prob- 
lems with these Bonnor-Ebert models. First, most star- 
less cores exhibit elongated shapes (see, e.g.. Fig. ^ and 
Fig. Ija) which are not consistent with spherical models 
beyond ~ 0.03 pc. Second, the inferred density contrasts 
(from center to edge) are generally larger (i.e., ^ 20-80 
- see Fig. ^d) than the maximum contrast of ~ 14 for a 
stable Bonnor-Ebert sphere (cf. § 1.1). Third, the effec- 
tive core temperature needed in these models (for ther- 
mal pressure gradients to balance self-gravity) is often 
significantly larger than both the average dust temper- 
atures measured with ISOPHOT (e.g. Ward-Thompson 
et al. 2002) and the gas temperatures measured in NH3 
(e.g. Lai et al., this volume). In the case of L1689B, for 
instance, the effective temperature of the Bonnor-Ebert 
fit shown in Fig. ^b is T^ff ~ 50 K, while the dust tem- 
perature observed with ISOPHOT is only Td ^ 11 K 
(Ward-Thompson et al. 2002). Finally, the physical pro- 
cess responsible for bounding the cores at some external 
pressure is unclear. These arguments suggest that prestel- 
lar cores cannot be described as simple isothermal hydro- 
static structures and are either already contracting (see 
§3.1 below) or experiencing extra support from static or 
turbulent magnetic fields (e.g. Curry & McKee 2000). 

One way of accounting for large density contrasts and 
high effective temperatures is to consider models of cores 
threaded by a static magnetic field and evolving through 
ambipolar diffusion (e.g. Ciolek & Mouschovias 1994, Basu 
& Mouschovias 1995). In these models, at any given time 
prior to protostar formation, the cores are expected to 
feature a uniform-density central region whose size cor- 
responds to the instantaneous Jeans length. This agrees 
well with the characteristics of the flat inner regions seen 
in starless cores. Furthermore, the observed sharp edges 
(e.g. Fig. ||b) are consistent with the model predictions 
(shortly) after the formation of a magnetically supercrit- 
ical core. Physically, this is because when a supercriti- 
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cal core forms, it collapses dynamically inward, while the 
outer, subcritical envelope is still efficiently supported by 
the magnetic field and remains essentially "held in place" . 
As a result, a steep density profile develops at the outer 
boundary of the supercritical core (see Fig. 8 of Basu & 
Mouschovias 1995). 

A problem, however, with these ambipolar diffusion mod- 
els involving only a static magnetic field is that they re- 
quire fairly large field strengths (~ 30-100 /^G, see Bac- 
mann et al. 2000), which seem to exceed the (few) exist- 
ing Zeeman measurements for low-mass dense cores (e.g. 
Crutcher 1999, Crutcher & Troland 2000). It is possible 
that more elaborate versions of the models, incorporating 
the effects of a non-static, turbulent magnetic field in the 
outer parts of the cores and in the ambient cloud, would 
be more satisfactory and could also account for the fil- 
amentary shapes often seen on large ( ^ 0.25 pc) scales 
(cf. Curry 2000, 2002 and Jones & Basu 2002). 

2.3. Temperature distribution in starless cores 

The models discussed in § 2.2 above assume that prestel- 
lar cores are isothermal, which is quite a good first ap- 
proximation (e.g. Larson 1969, Tohline 1982). In actual 
fact, however, there are good reasons to believe that the 
central regions of starless cores are somewhat cooler than 
their outer regions. Indeed, the results of ISO imaging 
in the mid-IR (Bacmann et al. 2000) and far-IR (Ward- 
Thompson et al. 2002) are consistent with the idea that 
starless cores are heated only from outside by the local 
interstellar radiation field (ISRF): the amount of energy 
radiated by such cores in the far-IR is roughly equal to 
the fraction of ISRF energy absorbed at shorter (e.g. mid- 
IR) wavelengths, with no evidence for any central heating 
source (Ward-Thompson et al. 2002). In such a situation, 
recent dust radiative transfer calculations (Evans et al. 
2001; Zucconi et al. 2001) predict that there should be 
a positive temperature gradient from the cores' centers 
(with Td as low as ^ 5-7 K) to their edges (at Td ^ 15 K). 
[Note that, in the dense (nnj > 3 x lO^cm"'^) inner re- 
gions, one expects the temperature of the gas to be well 
coupled to that of the dust (e.g. Do ty fc Neufeld 1997).] 
The first calculations published by Evans et al. (200T)| 
and Zucconi et al. (2001)| assumed rather a simplistic in- 
put radiation field, based on current estimates of the av- 
erage radiation field in the solar neighborhood (e.g. Black 
1994). We have more recently carried out similar calcula- 
tions in which we use the levels of the diffuse mid-IR and 
far-IR backgrounds observed toward the cores to make 
more realistic estimates of the effective radiation field di- 



rectly impinging on their surfaces. Fig. ^ gives an illus- 
tration for the prestellar core L1689B already discussed 
in § 2.1 (see FigJ|). In this case, the radiation field differs 
(in strength and spectrum) from the Black standard ra- 
diation field for several reasons: (1) due to the presence 
of early-type stars such as HD 147889 (B2 V) only a few 
parsecs away from the Ophiuchus cloud (e.g. Liseau et al. 
1999), the ambient far-UV field is significantly stronger 
(i.e.. Go ^ 12) than the solar neighborhood average; 
(2) this strong FUV field in turn excites small PAH-like 
grains in the outer [Ay ^ 1) layers of the cloud which 
are responsible for a high level of diffuse mid-IR emission 
distributed in several spectral bands (cf. Bernard et al. 
1993 and Fig. |a); (3) the L1689B core is itself embedded 
at some depth {Ay ~ 1) inside the cloud, so that a good 
fraction of the external FUV radiation is effectively repro- 
cessed to longer wavelengths. Consequently, the effective 
radiation field for L1689B is about one order of magni- 
tude stronger in both the mid-IR and the far-IR than 
the Black radiation field (cf. Fig. ^), in agreement with 
the strong mid-IR and far-IR backgrounds measured by 
Bacmann et al. (2000) and Ward-Thompson et al. (2002) 
with ISOCAM and ISOPHOT (- 4 MJy/sr at 7 ^im and 
- 120 MJy/sr at 200 fim, as opposed to - 0.25 MJy/sr 
and ~ 30 MJy/sr, respectively, in the Black radiation 
field). Based on this estimate of the effective radiation 
field, we have used the MODUST radiative transfer code 
(e.g. Bouwman 2001, Kemper et al. 2001) to compute the 
temperature distribution in the L1689B core for several 
plausible density distributions consistent with the 1.3 mm 
continuum emission map of Andre et al. (1996) and the 
7 fim absorption map of Bacmann et al. (2000). The re- 
sults are illustrated in Fig. |^. 

Our calculations confirm some of the conclusions of 



Evans et al. (2001) and Zucconi et al. (2001) , namely 



that the temperature reaches a minimum < 10 K in the 
centers of prestellar cores and that the central tempera- 
ture depends primarily on the central optical depth (di- 
rectly related to the degree of shielding from the external 
ISRF). However, we find that the minimum temperature 



may not be as low as reported by Evans et al. (2001) 



and Zucconi et al. (2001) and that the shape of the tem- 
perature profile may be quite different: it can be seen 
in Fig. Ijb that the major drop in temperature occurs 
in the outer parts of the core in our preferred model of 
L1689B (red curve), while it occurs closer to the center 
in the Evans et al. (2001) model (dashed curve). We also 
stress that all of the models we have run have a relatively 
high mass-averaged dust temperature of ~ 11 K, com- 
parable to the temperature estimated by fitting a grey- 
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Figures. Averaged 1.3 mm radial intensity profiles of (a) L1544 (prestellar) and L1527 (Class 0) in Taurus (d = 
140 pc), and (b) HH211-MM and L14-48-N (Class Os) in Perseus (d — 350 pc), compared to a synthetic profile, 
simulated for a singular isothermal sphere (SIS) model at T — 10 K (see Appendix of Motte & Andre 2001 for details 
about such simulations). The profile of the 30m beam is shown. An approximate column density scale, calculated 
assuming representative dust properties (Td — 15 K, K1.3 — 0.0075 cm^g^^J is also indicated. (Adapted from Ward- 
Thompson et al. 1999 and Motte & Andre 2001.) 



body to the global SED (cf. Fig. 13 of Ward- Thompson et 
al. 2002). (Note that Doty & Palotti 2002 have indepen- 
dently reached a similar conclusion.) Furthermore, if our 
assumptions about the effective radiation field and the 
dust opacity law (see Fig. ||b) are correct, then the radius 
of the flat inner region in the density profile of L1689B 
must lie between Rfiat — 3300 AU (smallest value con- 
sistent with the 1.3 mm constraints of Andre et al. 1996) 
and Rfiat = 6600 AU (value found by Bacmann et al. 
2000 assuming a central core temperature of 12.5 K - see 
§ 2.1). A value as small as Rfiat = 2000 AU (cf. Evans 
et al. 2001) is ruled out according to our analysis. Thus, 
conclusions derived from (sub)millimeter emission maps 
on the density structure of prestellar cores under the sim- 
plifying hypothesis that the dust temperature is uniform 
(e.g. Ward-Thompson et al. 1999) are probably not seri- 
ously flawed. 

In our opinion, the final word on the structure of 
prestellar cores will come from high-resolution mapping 
at far-IR and submillimeter wavelengths with future space- 
borne telescopes such as Herschel to be launched by ESA 
in 2007 (e.g. Pilbratt et al. 2001). Using Herschel im- 
ages at 75-500 /xm (in combination with ground-based 



dust continuum mapping at longer submillimeter wave- 
lengths) to construct SED maps for at least the nearest, 
spatially resolved prestellar cores, it will be possible to 
reconstruct their intrinsic temperature and density dis- 
tributions simultaneously as illustrated in Fig. ^. Such 
multi-band mapping observations will also set strong con- 
straints on the dust emissivity properties (e.g. the opacity 
index /3) and their spatial variation across dense cores. 

2.4. Structure of protostellar envelopes 

In the case of protostellar envelopes, the accretion lumi- 
nosity released close to the central protostar provides an 
additional source of heating compared to starless cores. In 
practice, however, the thermal structure of the envelopes 
surrounding low- luminosity protostars (Lbo; ~ 1-10 i©), 
including most Taurus embedded YSOs, is strongly af- 
fected by external heating from the ambient ISRF (Shirley 
et al. 2002), and the simple isothermal assumption (with 
T ~ 10 K) is often not too bad over a wide range of radii 
(r - 3000-15000 AU - cf. Motte & Andre 2001). 

The density structure of YSO envelopes has been stud- 
ied primarily with the (sub)millimeter emission approach 
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(although see Harvey et al. 2001 for a recent extinction 
study). In contrast to prestellar cores, protostellar en- 
velopes are always found to be strongly centrally con- 
densed (see, e.g.. Fig. ^d) and do not exhibit any marked 
inner flattening in their radial (column) density profiles. 
This is illustrated in Fig. |^ which shows the circularly- 
averaged radial intensity profiles measured at 1.3 mm 
for the circumstellar envelopes of the Class objects 
L1527 in Taurus and HH211-MM and L1448-N in Perseus 
(Motte & Andr 2001). (Note that the inner disk contribu- 
tion has been subtracted from the profiles based on high- 
resolution interferometric measurements with the IRAM 
Plateau de Bure interferometer.) Apart from a small over- 
density factor 2), the profile of the L1527 envelope 
closely follows the shape of the SIS model profile (see 
Fig. ||a), and is significantly steeper than the profile of 
the prestellar core LI 544 at small radii. 

More generally, several (sub)millimeter continuum stu- 
dies indicate that protostellar envelopes in regions of iso- 
lated star formation such as Taurus have radial density 
gradients consistent with p(r) cx r^^ with p 1.5-2 over 
more than ~ 10000-15000 AU in radius (e.g. Ladd et 
al. 1991, Chandler & Richer 2000, Hogerheijde & SandeU 
2000, Shirley et al. 2000, Motte & Andr 2001, J0rgensen 
et al. 2002). The density gradient estimated for (Class 
and Class I) protostars (i.e., p ^ 1.5-2) agrees with most 
collapse models which predict a value of p between 2 and 
1.5 during the protostellar accretion phase (before and 
after the passage of the collapse expansion wave, respec- 
tively - see § 1.1 above). Furthermore, the median en- 
velope mass measured by Motte & Andre (2001) for the 
bona-fide protostars of Taurus [Afe„„(r < A200 AU) ^ 
0.3 Mq] is roughly consistent with the predictions of the 
inside-out collapse model of Shu (1977) for ~ 10^ yr-old 
protostars. 

The situation is markedly different in regions where 
stars form in tight groups or clusters, such as Serpens, 
Perseus, and the p Oph main cloud. In this case, the 
observed envelopes are clearly not scale-free: they merge 
with dense cores, other envelopes, and/or the diffuse am- 
bient cloud at a finite radius R„^t ^ 5 000 AU (Motte, 
Andre, k Neri 1998; Motte & Andre 2001; see also Fig. 5 
of Mundy, Looney, & Welch 2000). Moreover, as can be 
seen in Fig. |b for HH211-MM and L1448-N, the Class 
envelopes mapped in cluster-forming regions arc found 
to be 3 to 12 times denser than the SIS at T= 10 K 
(which is the typical gas kinetic temperature expected in 
these clouds prior to massive star formation - see Gold- 
smith & Langer 1978, Evans 1999, Goldsmith 2001, and 
Fig. ||b) . Static magnetic fields can only account for mod- 



erate overdensity factors ( ^ 2) compared to an unmag- 
netized SIS (see Li & Shu 1997). Turbulence could in 
principle contribute to the support of the initial dense 
core (e.g. Myers & Fuller 1992, Mardones et al. 1997). 
However, the small-scale condensations corresponding to 
the precursors of protostars in star-forming clusters (e.g. 
Motte et al. 1998) appear to be largely devoid of tur- 
bulence (see § 3.1 below). A more likely explanation is 
that, in such regions, protostellar collapse is induced by 
strong external disturbances and starts from non-singular 
initial conditions, resulting in a non- equilibrium density 
configuration with a large overdensity factor ( ^ 10) and 
significant inward velocities near t ~ (see also § 4 and 
Fig. P below). 

3. Velocity Structure 

3.1. Internal motions in starless cores 

On large ( ^ pc) scales, the spectral line profiles/widths 
observed in molecular clouds are indicative of supersonic 
turbulent motions (e.g. Falgarone & Phillips 1990). The 
nonthermal component of the line width, catt, roughly 
follows the scaling law gnt oc (e.g. Fuller & My- 

ers 1992), which is believed to arise from MHD turbu- 
lence (e.g. Mouschovias & Psaltis 1995). On intermediate 
( ^ 0.1 pc) scales, there are strong differences between 
the "low-mass" cores observed in regions of distributed 
star formation (e.g. Taurus), whose line widths are dom- 
inated by thermal motions (e.g. Benson & Myers 1989), 
and the "massive", cluster-forming cores of, e.g., Ophi- 
uchus and Orion, whose line widths are still dominated 
by turbulent motions (e.g. Caselli & Myers 1995). On 
small ( ^ 0.03 pc) scales, however, the prestellar conden- 
sations of the Ophiuchus, Serpens, Perseus, and Orion 
cluster-forming regions are always characterized by very 
narrow line widths (e.g. Belloche et al. 2001, Myers 2001), 
which is reminiscent of the thermal cores of Taurus. In the 
p Oph protocluster, for instance, the nonthermal velocity 
dispersion is about half the thermal velocity dispersion of 
H2 {(Tnt/ctt ~ 0.7) toward the starless condensations of 
the dense cores Oph Bl, C, E, F (Belloche et al. 2001). 
This indicates that, even in cluster-forming clouds, the 
initial conditions for individual protostellar collapse are 
"coherent" (cf. Goodman et al. 1998) and essentially free 
of turbulence. The dissipation, on small (< 0.1 pc) scales, 
of a significant fraction of the turbulent motions observed 
on large scales is thus a prerequesite for the formation of 
prestellar condensations (e.g. Nakano 1998, Myers 1998, 
Falgarone & Pety 2001). Apart from low levels of turbu- 
lence, low-mass prestellar cores are also characterized by 
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small rotation rates in general {Q, ~ 10"^'* — 10~^^ rad s~^ 
- e.g. Goodman et al. 1993), and subsonic, extended in- 
fall motions (e.g. Lee, Myers, & Tafalla 2001 - see also 
Evans, this volume). 

3.2. Rotation and infall in Class objects 

The velocity structure of only a few Class objects has 
been studied in detail up to now. Here, we mostly discuss 
the example of the nearby (d = 140 pc) and relatively iso- 
lated Class object IRAM 04191+1522 (IRAM 04191 for 
short) found by Andre, Motte, & Bacmann (1999) in the 
vicinity of the Taurus Class I source IRAS 04191+1523 
(e.g. Tamura et al. 1991 - see Fig. |b). 

IRAM 04191 features a very high envelope mass to 
luminosity ratio {Mf^^°°^^ / Lboi ~ 3 Mq/Lq) and a low 
bolometric temperature (Tboi ~ 18 K), suggesting an age 
t 3 X 10"*^ yr since the beginning of the accretion phase 
(see Andre et al. 1999). It is also remarkable in that it has 
a very low bolometric luminosity {Lboi ~ 0.15 L©), and 
no or only a tiny accretion disk {Rdisk < 15 AU, Mdisk < 
2 X 10"^ Mq). Accordingly, IRAM 04191 is probably the 
youngest accreting protostar currently known in Taurus. 

IRAM 04191 is associated with a flattened circumstel- 
lar envelope, seen in the dust continuum (Fig. |l]b) and 
in dense gas tracers such as N2H+, C3II2, H^'^CO^, and 
DCO+ (Belloche et al. 2002). All of the maps taken in 
small optical depth lines show a clear velocity gradient 
across this envelope of ^ 7kms"^pc"^ increasing from 



west to east along the long axis, i.e., roughly perpendicu- 
lar to the outflow axis (see Fig. |^b and Fig. ||). Strongly 
suggestive of rotation, this velocity gradient corresponds 
to a mean angular velocity ~ 3 x 10"^'^ rad s~^ in 
the inner ~ 3000 AU region (after deprojection from a 
viewing angle i = 50°). Furthermore, this rotation does 
not occur in a rigid-body, but differential, fashion: the 
inner ~ 3500 AU-radius region seems to be rotating sig- 
nificantly faster than the outer parts of the core (see 
Fig. H). The "S" shape of the position- velocity diagram 
shown in Fig. ^ can be fitted with a rotation curve of the 
form Vrot{r) oc r°-^±OA ^ ^ 35QQ j^jj ^^^^ Vrot{r) oc 
j,-i.5±o.5 35QQ < r < 11000 AU, corresponding to 
il(r) c>c 7--o.9±o.4 g^^^ Q{r) oc 7--2.5±o.5^ respectively (see 
Fig. |a and Belloche et al. 2002 for details). 

Direct evidence for infall motions over a large portion 
of the IRAM 04191 envelope is also observed in optically 
thick lines such as CS(2-1), CS(3-2), H2CO(2i2 - In), 
and H2CO(3i2 — 2ii). These lines are double-peaked and 
skewed to the blue up to ^ 40" from source center (see 
Fig. 1^), which is indicative of infall motions up to a radius 
Rinf ~ 5000 AU from center (cf. Evans 1999 and Evans, 
this volume). Radiative transfer modeling confirms this 
view, suggesting a relatively flat infall velocity profile 
{Vinf ~ 0.1kms~i) for 3000 <r< 11000 AU and larger 
infall velocities scaling as Vmf oc r^°-^ for r ^ 3000 AU 
(see Fig. |b and Belloche et al. 2002 for details). The 
mass infall rate is estimated to be Mmf ^ 2 — 3xal/G^ 
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3 X 10"^ Mq yr-i (with - 0.15 - 0.2 km s^^ for T - 
6 — 10 K), roughly independent of radius. 

Another Class object whose kinematics has been 
quantified in detail recently is the IRAS 4 system in the 
NGC 1333 protocluster (Di Francesco et al. 2001). IRAS 4 
is a triple system (A/B/C) on scales of ^ 10000 AU, 
whose primary component (A) itself breaks up into a ~ 
600 AU binary (e.g. Looney et al. 2000). Using the IRAM 
Plateau de Bure interferometer, Di Francesco et al. (2001) 
observed inverse P-Cygni profiles in H2CO(3i2 — 2ii) to- 
ward IRAS 4 A and IRAS 4B, from which they derived 
very large mass infall rates of ~ 1.1 x 10"'' Mq yr"'^ and 
~ 3.7 X lO-'^ Mq yr-i at r - 2000 AU for A and B, 
respectively. Even if a warmer initial gas temperature (~ 
20 K) than in IRAM 04191 and some initial level of tur- 
bulence are accounted for (see Di Francesco et al. 2001), 
these values of Minf correspond to more than ^15 and 
~ 6 times the canonical a^f^/G value, respectively (where 
Ues ^ 0.3 kms~^ is the effective sound speed). These very 
high infall rates result both from very dense envelopes 
(see Motte & Andre 2001) and large, supersonic infall ve- 
locities (~ 0.68 km s^^ and ~ 0.47 km s"'^ for A and B, 
respectively, according to Di Francesco et al. 2001). Ev- 
idence of fast rotation, producing a velocity gradient as 
high as ^ 40 km s~^ pc"'^ , was also found by Di Francesco 
et al. (2001) toward IRAS 4A. 

4. Conclusions: Comparison with collapse models 

In the case of low-mass, isolated dense cores, the SIS 
model of Shu (1977) describes global features of the col- 
lapse (e.g. the mass infall rate) reasonably well and thus 
remains a useful, approximate guide. In detail, however, 
the extended infall velocity profiles observed in prestellar 
cores (e.g. Lee et al. 2001) and in the very young Class 
object IRAM 04191 (§ 3.2 above) are inconsistent with 
the pure inside-out collapse picture of Shu (1977). The 
shape of the density profiles observed in prestellar cores 
are well fitted by purely thermal Bonnor-Ebert sphere 
models, but the absolute values of the densities are sug- 
gestive of some additional magnetic support (§ 2.2). The 
observed infall velocities are also marginally consistent 
with isothermal collapse models starting from Bonnor- 
Ebert spheres (e.g. Foster & Chevalier 1993, Hennebelle 
et al. 2002), as such models tend to produce somewhat 
faster velocities. This suggests that the collapse of 'iso- 
lated' cores is essentially spontaneous and somehow mod- 
erated by magnetic effects in magnetized, probably not 
strictly isothermal Bonnor-Ebert cloudlets. Furthermore, 
the contrast seen in Fig. p| between a steeply declining 



rotation velocity profile and a flat infall velocity profile 
beyond ^ 3500 AU suggests that angular momentum is 
not conserved in the outer envelope of IRAM 04191. Such 
a behavior is very difficult to explain in the context of 
non-magnetic collapse models. In the presence of a rel- 
atively strong (~ 60 fiG) magnetic field, on the other 
hand, the outer envelope can be coupled to, and spun 
down by, the (large moment of inertia of the) ambient 
cloud (e.g. Basu & Mouschovias 1994). Based on a quali- 
tative comparison with the ambipolar diffusion models of 
Basu & Mouschovias, Belloche et al. (2002) propose that 
the rapidly rotating inner envelope of IRAM 04191 cor- 
responds to a magnetically supercritical core decoupling 
from an environment still supported by magnetic fields 
and strongly affected by magnetic braking. (Quantita- 
tively, however, the models published so far are rotating 
too slowly to fit the observations.) In this view, the inner 
^ 3500 AU radius envelope of IRAM 04191 would cor- 
respond to the effective mass reservoir (~ 0.5 Mq) from 
which the central star is being built. Moreover, compar- 
ison of these results with the rotational characteristics 
of other objects in Taurus (Ohashi et al. 1997) suggests 
that IRAM 04191 behaves in a typical manner and is sim- 
ply observed particularly soon after point mass formation 
(i.e., at t ^ 0). If this is correct, the masses of the stars 
forming in clouds such as Taurus may be largely deter- 
mined by magnetic decoupling effects. 

In protoclusters, by contrast, the large overdensity 
factors measured for Class envelopes compared to hy- 
drostatic isothermal structures (cf. § 2.4), as well as the 
fast supersonic infall velocities and very large infall rates 
observed in some cases (e.g. § 3.2), are inconsistent with 
self-initiated forms of collapse and require a strong ex- 
ternal influence. This point is illustrated in Fig. ^ with 
the results of recent SPH simulations by Hennebelle et 
al. (2002). These simulations follow the evolution of a 
Bonnor-Ebert sphere whose collapse has been induced by 
an increase in external pressure Pext- Large overdensity 
factors (compared to a SIS), together with supersonic 
infall velocities, and large infall rates ( ^ \0a^^ /G) are 
found near t = Q when (and only when) the increase 
in Pea;t IS strong and very rapid (e.g. Fig. ||b), result- 
ing in a violent compression wave. Such a violent col- 
lapse may be conducive to the formation of both mas- 
sive stars (through higher accretion rates) and multiple 
systems (when realistic, non-isotropic compressions are 
considered). Future high-resolution studies with the next 
generation of (sub)millimeter instruments (e.g., ALMA) 
will greatly help test this view and shed further light on 
the physics of collapse in cluster-forming regions. 
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Figures, (a) Effective radiation field impinging on 
L1689B (red, solid curve), as estimated from the observed 
levels of the mid-IR and far-IR backgrounds (blue tri- 
angles, see text). This is much stronger than the Black 
(1994) radiation field (black, dashed curve), (b) Model 
dust temperature profiles calculated for L1689B with 
MODUST under two sets of assumptions: the red, solid 
curve assumes the effective radiation field shown in (a), 
the density profile derived by Bacmann et al. (2000) (with 
Rfiat = 6600 AU), and the dust opacity law of Preibisch 
et al. (1993) with thin ice mantles in the core and no 
ice mantles in the ambient cloud; the black, dashed curve 
corresponds to the model favored by Evans et al. ( 2001 ) 
which has Rfiat — 2000 AU (for d = 160 pc), no ambient 
cloud, and assumes the Black (1994) radiation field xO.5 
and the dust opacity law of Ossenkopf & Henning (1994) 
with coagulation (0H5). 
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Figure J^. Simulated Herschel observations of the radial 
structure of the prestellar core L1689B. The top panel (a) 
shows two model intensity profiles, consistent with present 
millimeter data, as they could be observed by Herschel at 
250 iim. The upper (red) curve corresponds to the model 
we favor here based on our best estimate of the effec- 
tive radiation field (see Fig. while the lower curve is 
for the model proposed by Evans et al. ( 2001 ) based on 
the Black (1994) ISRF xO.5. Note the large differences 
in shape and absolute scaling between the two models. 
The bottom panel (b) shows the dust temperature distri- 
butions (solid curves) calculated with MODUST for these 
two models. The dashed curves (almost indistinguishable 
from the solid curves beyond 20" ) show the reconstructed 
temperature profiles that could be derived from Herschel 
maps in 6 bands between 75 /im and 500 /im (plus ground- 
based data at 850 fj,m and 1.3 mm), assuming the 3-D 
core geometry is known. The reconstruction becomes un- 
certain below the Herschel resolution limit (9" HP beam 
radius at 250 fim). 
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Figure 6. Position- velocity diagram along the major axis 
of the IRAM 04191 envelope (i.e., perpendicular to the 
flow) based on a C*^^S'(2 — 1) map taken at the IRAM 30m 
telescope (Belloche et al. 2002). Contours: 0.2 to 0.8 by 
0.2 K. The dots with error bars mark the observed velocity 
centroids. The solid curve shows the profile of a model 
with differential rotation (Vrot{r > 3500 AU) cx r"^'^^ 
also used in Fig. (see text). 
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shaded areas show the estimated domains where the mod- 
els match the CS and C^'^S observations reasonably well. 
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Figure 9. Density profiles ( solid curves ) obtained slightly 
after point mass formation (t ^ 10* yr) in SPH nu- 
merical simulations of the collapse of a nearly critical 
[{pc/ Pout)init ~ 12/ isothermal (T ~ 10 K) Bonnor- 
Ebert sphere induced by external compression (Hennebelle 
et al. 2002). For comparison, the dotted line shows the 
p cx r"^ profile of a SIS at 10 K. In (a), the collapse 
was initiated quasi- statically (by very slow compression 
with Pext/Pext = 20 X the initial sound crossing time 
Rinit/o..,) Obnd the density profile at t ^ is very similar 
to that of the SIS. By contrast, in (b), the collapse was in- 
duced by a very rapid increase in external pressure (with 
Pext/Pext — 0.03 X Rinit/cLs), resulting in much larger 
densities around t ^ 0. 



